Self-phase modulation of a single-cycle terahertz pulse by nonlinear free-carrier 

response in a semiconductor 
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We demonstrate the self-phase modulation (SPM) of a single-cycle THz pulse in a semiconductor, 
using bulk n-GaAs as a model system. The SPM arises from the heating of free electrons in the 
electric field of the THz pulse, leading to an ultrafast reduction of the plasma frequency, and hence 
to a strong modification of the THz-range dielectric function of the material. THz SPM is observed 
directly in the time domain. In the frequency domain it corresponds to a strong frequency-dependent 
refractive index nonlinearity of n-GaAs, found to be both positive and negative within the broad 
THz pulse spectrum, with the zero-crossing point defined by the electron momentum relaxation 
rate. We also observed the nonlinear spectral broadening and compression of the THz pulse. 



Since its first demonstration [l[ more than five decades 
ago, nonlinear optics (NLO) has revolutionized the field 
of photonics. Until recently, ultrafast NLO was largely 
confined to the visible and infrared spectral ranges. How- 
ever, with the emergence of sources generating high- 
energy ultrafast single-cycle terahertz (THz) pulses [2- 
IH, more and more NLO studies are being conducted in 
the THz range. Nonlinear dynamics in semiconductors 
in strong THz fields is a particularly rich area of research 
due to good coupling between the THz fields and free 
carriers, and due to the possibility of inducing strong po- 
larization effects in semiconductors with strong THz sig- 
nals . Many interesting ultrafast phenomena have been 
investigated recently, such as interaction of strong THz 
fields with excitons [g, 0] , THz-induced impact ionization 
of electrons in narrow-gap semiconductors [8|, Q , ultrafast 
THz-induced electron dynamics in complex conduction 
bands [10 13], and femtosecond all-optical switching in 
nano-scale structures with THz signals [13, [HI • 

The majority of studies in the present-day ultrafast 
THz NLO are focused on achieving a better understand- 
ing and control of nonlinear processes in (mostly semi- 
conductor) materials using THz light. There are, how- 
ever, very few reports aiming at control and manipulation 
of THz light making use of the material's nonlinear re- 
sponse. Yet, the ultrafast THz signals are usually single- 
cycle waveforms, which can be experimentally sampled 
with a time resolution much better than one optical cycle, 
providing the complete time-domain information about 
the optical signal - its amplitude and phase 
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we will show in this Letter, this makes the THz pulses 
ideal tools for the direct observation of both i) the THz- 
specific NLO phenomena, and it) the general NLO phe- 
nomena occurring in the single-cycle optical regime. 

Single-cycle pulses can not be described by the tra- 
ditional carrier-envelope approximation (l8j |. Irrespec- 
tive of the frequency range to which they belong, they 
inherently have an extremely broad spectral bandwidth 
covering many octaves of frequencies. Hence, the propa- 
gation of such signals in the medium is governed by the 



complex- valued material dielectric function s(uj) cover- 
ing a very broad spectral range. An optical nonlinearity 
- the near-instantaneous modification of this dielectric 
function due to a nonlinear light-matter interaction oc- 
curring in the single-cycle regime, corresponds to a mod- 
ification of the optical properties of matter over an ex- 
tremely broad range of frequencies. We find the effect 
of such a sudden change of dielectric environment on the 
propagating single-cycle waveform to be quite intriguing, 
which is the motivation behind the present work. 

In this Letter we report on the study of a fundamental 
nonlinear optical effect - the self phase modulation (SPM) 
[H S3, of a single-cycle THz pulse. The SPM occurs 
due to the change of the (real part of) refractive index 
n(cj) = Re y/£(u}) of the material as a result of light- 
matter interaction. As an efficient nonlinear medium for 
the THz range we used an n-doped bulk semiconduc- 
tor GaAs, whose optical nonlinearity arises from the re- 
sponse of free electrons to the THz fields. The n-GaAs 
was chosen because it is a well known, model system; 
however, any other doped or photoexcited semiconduc- 
tor with a complex band structure will also exhibit the 
THz-range SPM according to the mechanism described in 
this Letter. The SPM was observed by us directly in the 
time domain using THz time-domain spectroscopy (TDS) 
HI, l3], and was analyzed in the frequency domain. In 



particular, we found that the nonlinear refractive index 
An(w) of n-GaAs is strongly frequency-dependent: it ex- 
hibits both positive and negative nonlinearity within the 
broad spectral bandwidth of a single-cycle THz pulse, 
with the sign of nonlinearity determined by the the rate 
of electron momentum relaxation in the semiconductor. 
We also observed the nonlinear spectral broadening and 
compression of the THz pulses in the material. 

In this Letter we show that the THz-range optical non- 
linearity of a doped or photoexcited semiconductor is 
based on the modification of its plasma frequency by the 
carrier heating in the THz field. Semiconductors with 
free carriers are good absorbers of THz radiation, and 
their THz-range complex- valued dielectric function i(u>) 



2 







m x > m r \ 




\ / m L > m r 




n - GaAs ™ rn r 




high THz absorption 


low THz absorption 




0.5 t 



10 12 14 



FIG. 1. (a) Mechanism of THz nonlinearity in n-GaAs, based 
on carrier heating in the electric field of a THz pulse, lead- 
ing to a decrease of plasma frequency ui p , and saturation of 
THz absorption, (b) Temporal dynamics of THz nonlinearity 
in n-GaAs, observed in a transmission-mode TPTP measure- 
ment. Near-instantaneous bleaching of the sample via carrier- 
heating in the THz field is followed by a slower recovery due 
to electron cooling. 



is well described by the basic Drude plasma model [2 iL |22 1 

£(uj) = (n + iac/2ui) 2 = Edc — ^p/(<^ 2 — iuj/r r ) (1) 



or by one of its extensions [23, [24{ . Here n and a are the 
frequency-dependent refractive index and power absorp- 
tion coefficient, respectively, r r is a carrier momentum 
relaxation time, and e ( ic is the background "static" di- 
electric constant of the semiconductor in the absence of 
free carriers 25|. uj p = (Ne 2 /eom) 1 / 2 is the plasma fre- 



quency, where e is the elementary charge, N is the free 
carrier density, eo is the vacuum permittivity, and m is 
the carrier effective mass. 

It directly follows from Eq. (TTJ) that the change in 
plasma frequency uj p resulting from the THz light - mat- 
ter interaction will lead to a change in the THz absorption 
a and the refractive index n of the semiconductor, i.e. to 
nonlinear absorption and SPM. In the limiting case of 
ujp — >• the THz absorption will vanish, and the THz re- 
fractive index will be defined by the background dielectric 
constant n(cu p — > 0) = yfEdc- The only mechanism to re- 
duce the plasma frequency in a doped semiconductor (i.e. 
where N is fixed), is by increasing the carrier effective 
mass to. In n-doped materials (e.g. n-GaAs) this can be 
achieved by electron heating in a sufficiently strong elec- 
tric field - by redistribution of the electron population 
in energy-momentum space of the conduction band from 
low-momentum, small-effective-mass states around the 
bottom of the T- valley to high-momentum, large-effective 
mass states in the satellite valleys and to the states of 
strong nonparabolicity within the T-valley (see e.g. j2(| 
and references therein). In p-doped materials the hole 
scattering between the hole subbands can be utilized. 

Recently the carrier heating in the ponderomotive po- 
tential created by the electric field of a single-cycle THz 
pulse was demonstrated. In [Io| it was shown 
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FIG. 2. THz SPM in n-GaAs in the time domain. Peak 
Ethz of reference THz pulse in (a) is 292 kV/cm, and in 
(c) is 9 kV/cm. (b) Reference and sample pulses normalized 
to their maxima, for the whole range of peak Ethz of 9 - 
292 kV/cm. Black dashed line: estimated group delay of 
the sample pulse. (d,e) Normalized reference (d) and sample 
(e) pulses for the cases of peak Ethz of 9 kV/cm and 292 
kV/cm, demonstrating the SPM in the sample, (f) Amplitude 
frequency-integrated transmission coefficient, and estimated 
group delay of a reference pulse as a function of peak Ethz- 



that right after the excitation of n-GaAs by a THz pulse 
with peak electric field Ethz — 100 kV/cm, the effec- 
tive mass of electrons initially occupying the bottom 
of T-valley (m-p = 0.063mo) drastically increases and 
reaches the values of m ~ 0.3mo, more typical for the 
L— and X~ valley states in GaAs (mj, = 0.23mo, and 



mx = 0.43mo) [too is the electron rest mass] |26|, |27 |. 
The THz saturable absorption based on this mechanism, 
shown in Fig. [IJa), was observed by us in various n- 
doped semiconductors (28|. The temporal response of 
THz nonlinearity of n-GaAs is near-instantaneous, as 
shown in Fig. [TJb). This measurement is a result of 
THz pump - THz probe (TPTP) spectroscopy [13, El El 
where our n-GaAs sample with the thickness d = 0.4 mm, 
Af = 5x 10 15 cm~ 3 , and r r = 94 fs was excited by a THz 
pump pulse with peak Ethz — 200 kV/cm. The near- 
instantaneous bleaching of the sample is followed by a 
slower recovery as the electrons return to the bottom of 
the T- valley [291 ] - This, and all other measurements in 
this work were performed at room temperature. 

We observed the SPM of single-cycle THz pulses in 
our n-GaAs sample in a nonlinear THz TDS experiment 
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FIG. 3. THz SPM in n-GaAs in spectral domain. Nonlin- 
ear contribution to refractive index An, and corresponding 
nonlinear phase delay as a function of frequency and peak 
Ethz in the range 9 - 292 kV/cm. Blue arrows show the 
positive and negative refractive index nonlinearity with in- 
creasing peak Ethz- The amplitude spectrum of reference 
THz pulse is shown as a background. 



set up in transmission configuration [28j . The frequency 
spectrum of the THz pulses, generated by optimized opti- 
cal rectification of 800 nm, 100 fs lasers pulses in lithium 
niobate (LN) [2|-|4|, covered the range 0.1-3 THz, thus 
containing many octaves of frequencies. The peak elec- 
tric field Ethz in the generated THz pulses was adjusted 
in the range of 9 - 292 kV/cm using a pair of wire- 
grid polarizerspositioned after the THz emitter in the 
spectrometer [2|, |28j . THz transmission through the po- 
larizers is nominally frequency-independent. The THz 
pulses were recorded in the time domain using free-space 
electro-optic sampling 1(| 17 1. 

As shown in Fig. [21 for each value of peak Ethz in the 
generated THz pulse a reference (travelling through the 
optical path in the spectrometer in vacuum) and a sample 
(travelling through the same geometrical path, but with 
the sample inserted into it) THz pulses [lg, [l7|, centered 
around the time delays of ps and 3.3 ps, respectively, 
were recorded. A clear THz saturable absorption (i.e. 
increase in THz transmission) [28| in our n-GaAs sample 
with growth in peak electric field Ethz is observed. In 
Figs. [U[a) and (c) the measured reference and sample 
THz pulses for the two extreme values of peak Ethz, 
292 kV/cm and 9 kV/cm, are shown as examples of raw 
data. Here and below all mentioned peak values of Ethz 
correspond to that of the reference THz pulse, and are 
reported as measured at the position of the sample. 

In Fig. [H[b) a two-dimensional plot shows the refer- 
ence and sample pulses normalized to their respective 
peak values, for the whole range of peak Ethz in our 
experiments. As apparent from this plot, the time de- 
lay between the reference and the sample THz pulses 
in n-GaAs increases with increasing peak Ethz, which 
is a direct signature of SPM. The dashed black line in 
Fig. [U(b) shows the estimated group delay (i.e." delay 
as a whole") [28| of a THz pulse in the material as a 
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FIG. 4. Nonlinear spectral broadening and compression of a 
THz pulse in n-GaAs. (a) The ratio of effective bandwidths of 
the sample and reference pulses, as a function of peak Ethz ■ 
(b-d) Examples of spectral dynamics: normalized amplitude 
spectra of reference and sample THz pulses measured at peak 
THz field values of of 30 kV/cm, 145 kV/cm, and 198 kV/cm, 
indicated by arrows in (a). 



function of peak Ethz- In Figs. [2{d, e) the normalized 
reference and sample pulses for the two extreme values 
of peak Ethz of 9 kV/cm and 292 kV/cm are shown. 
The difference between the two reference waveforms in 
Fig. [5] (d) is fairly small: there is no observable time 
shift, and there is only a very small change in the pulse 
shapes, which is due to non-ideality of the THz wire-grid 
polarizers used for adjustment of the peak THz field. The 
difference between the two sample waveforms in Fig. [5] 
(e) is, however, quite noticeable. The sample waveform 
corresponding to a stronger Ethz clearly experiences a 
larger time delay, as well as a certain reshaping. Fig. [2] 
is therefore the direct time-domain observation of THz 
SPM occurring in a single-cycle regime. In Fig. [2] (f) 
we show the connection between the saturable absorp- 
tion and SPM in our sample, by plotting both the THz 
frequency-integrated amplitude transmission coefficient 
and an estimated group delay of the sample pulse as a 
function of peak Ethz- The observed increase in the 
transmission coefficient is more than two-fold, and the 
increase in the group delay exceeds 150 fs, making it a 
noticeable fraction of the duration of the THz cycle. 

We now transfer to the frequency domain by using the 
Fourier transforms of the measured THz pulses from Fig. 
[H 16l Il7|. In Fig. [3] we show the SPM as a nonlinear, 



THz-field dependent contribution to the refractive index 
of n-GaAs, calculated as a correction to the refractive 
index spectrum measured at the weakest peak Ethz , 
An(w) = n(w, Ethz) — n(u),9 kV/cm). A very strong, 
and also frequency-dependent, dynamics of An(w) is ob- 
served with increase in peak Ethz from 9 kV/cm to 292 
kV/cm. Interestingly, the refractive index nonlinearity is 
found to be both negative and positive within the band- 
width of the same single-cycle THz pulse, making it quite 
a unique situation in nonlinear optics (though quite un- 
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FIG. 5. A frame from the movie [33| showing the evolution 
of the reference and sample THz pulses in time domain, and 
their normalized Wigner-Ville spectrograms, with increase in 
peak Ethz- 



derstandable given the ultra-broadband nature of single- 
cycle pulses). These negative and positive trends in the 
dynamics of An at the frequencies below and above ap- 
proximately 0.5 THz, with increase in peak Ethz, are 
shown with blue arrows in Fig. [3] The maximum refrac- 
tive index changes in our measurements were as large 
as An — -0.13 — +0.08, in the spectral range 0.3 - 2.4 
THz. This corresponds to the nonlinear phase delay 
At-0 = And/c in the range of -173 — K107 fs acquired by 
the THz signal at different frequencies in our 0.4 - mm 
thick sample, also shown in Fig. |3] 

In our specific case most of the spectrum of the THz 
pulse (shown as grey area in Fig. [3]) was in the spectral 
range of positive index nonlinearity, which results in the 
increasing group delay of the THz pulse at higher peak 
Ethz (see Fig. [2]). However, if the THz pulse spectrum 
belonged to the range below 0.5 THz, the decrease in 
group delay with increase in Ethz would occur instead. 

Inherent to SPM is spectral reshaping of the optical 
signal 35, 3(|, which can lead to both spectral broaden- 
ing and compression within the same NLO system (see 
e.g. [I(J). In Fig. |4|a) it is illustrated by the dependency 
of the ratio of effective spectral bandwidths of the sample 



and reference pulses Aiu sam / Aui re f on peak Ethz [31 1. 
If this ratio is greater (smaller) than unity, then relative 
spectral broadening (compression) of a THz pulse in a 
sample takes place. We observed a clear spectral breath- 
ing by approx. -3 — +7 % of a THz pulse in our sam- 
ple. In Figs. |U(b-d) the normalized amplitude spectra 
of reference and sample THz pulses at selected values of 
peak Ethz are shown, exemplifying the spectral broad- 
ening, compression, and virtually no change, respectively, 
of THz pulses in n-GaAs. In general, the observed spec- 
tral breathing will depend on the combination of pulse 
chirp, index nonlinearity, and frequency-dependent non- 
linear absorption, which can only be quantified using ad- 
vanced numerical methods 0, [l3[ ■ 

A cognitive correlation between the THz pulse trans- 
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FIG. 6. Blue arrows: trends in (a) refractive index and (b) 
absorption THz spectra of n-GaAs with decrease in uj p . Spec- 
tra measured in linear regime (symbols), and calculated us- 
ing a DS model with iV = 5 x 10 15 cm" 3 , r r = 94 fs, and 
cos = —0.15 for the electron effective mass of Y- (red solid 
line), L- (blue dashed line), and X— (green dotted line) val- 
leys. Amplitude spectrum of a reference THz pulse is shown in 
the background of (a). Black arrows indicate the correspond- 
ing u) p values. Horizontal line in (a) corresponds to uj p — 0, 
when n = Jedc = 3.6. 



mission and delay as a function of peak Ethz, also de- 



picting the evolution of THz pulse spectrograms [32j, is 
visualized by a movie, which can be found under the Ref. 
33] . An example of the movie frame is shown in Fig. [5] 

We note that the link between the refractive index and 
absorption (or gain) nonlinearities is fundamental, and is 
used e.g. in semiconductor slow/fast light devices for the 
telecom range (see e.g. [34], HH). We also note that the 
reshaping of a THz signal generated in a LN crystal under 
variable optical pump strength was reported and assigned 
to THz SPM in Q, but later studies clearly related this 
effect to the optical-range nonlinearity of the laser pump 
pulse in LN [37J • The reshaping by SPM of a THz pulse 
in InSb undergoing impact ionization was reported 
however, without the systematic analysis. 

Now we illustrate the observed trends in nonlinear ab- 
sorption and refractive index by simple modeling. In Fig. 
|6] the measured linear THz spectra of our sample are 
shown (symbols) along with the Drude fit using the pa- 
rameters N and r r described above, n = ^/Sdc = 3.6 
[38| . and m e — mr (solid lines). For this fit we used a 
Drude-Smith (DS) model incorporating the effective car- 
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rier localization into an ideal Drude model 
The small localization parameter cds — —0.15 used here 
shows that our sample is very close to an ideal Drude 
plasma system described by Eq. [TJ In the same figure we 
show the calculated n and a spectra for the same sam- 
ple parameters as before, but using the heavier effective 
masses m e = mj, and m e = mx of the satellite valleys 
[27| | (dashed and dotted lines, respectively). Correspond- 
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ing values for oj p are shown in the figure using black ar- 
rows. It is clear that with decrease in uj p the absorption 
decreases monotonously. The trend in refractive index 
with decreasing u p is the decrease from larger values at 
lower frequencies (< 0.5 THz), and increase from lower 
values at higher frequencies, approaching the background 
value of n = y/£dc at o; p = 0. These trends in the dy- 
namics of a and n, indicated by the blue arrows, exactly 
match the observed saturable absorption (Fig. H{f)), and 
the behavior of refractive index nonlinearity (Fig. [3J , in- 
cluding its opposite signs at lower and higher frequencies. 
The large value of refractive index at lower frequencies in 
a metal or doped semiconductor directly follows from Eq. 
[TJ it is known as Hagen- Rubens regime,and is character- 
istic for the frequency range ui -C r^ 1 [22[. The spectral 



position of zero index nonlinearity is thus defined by the 
electron momentum relaxation rate r" 1 , and can be the- 
oretically predicted with reasonable accuracy, as shown 
by Figs. E] and EJa). 

In conclusion, we have studied the SPM of a single- 
cycle THz pulse in a semiconductor, whose THz-range 
optical nonlinearity is based on the carrier heating in 
THz fields, and is adequately described by the Drude 
plasma model. 
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